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Abstract 

We calculate the differential cross section for real photon electroproduction off spinless hadron 
which Sevres as a main probe of the hadrons structure via the concept of generalized parton 
distributions. Compared to previously available computations performed with twist-three power 
accuracy, we exactly accounted for all kinematical effects in hadron mass and momentum transfer 
which arise from leptonic helicity amplitudes. We performed numerical studies of these kinemat- 
ical effects and demonstrated that in the valence quark region and rather low virtualities of the 
hard photon which sets the factorization scale, the available approximate results significantly 
overestimate the cross section rates in comparison to exact formulas. 



1 Introduction 

The qualitative mechanism which binds hadronic constituents together in a bound state is well 
understood within the framework of QCD via the formation of coUimated gluon flux tubes be- 
tween quarks. This picture is confirmed by lattice gauge theory simulations which aim at a 
quantitative exploration of the hadronic structure. However, physical observables which can 
be measured in high-energy experiments necessarily involve correlation functions of elementary 
quark and gluon fields separated by light-like distances and thus evade straightforward use of 
euclidean lattice tools. Therefore, in the lack of other analytical/numerical first-principle tech- 
niques, phenomenological analyses of experimental data are currently the only viable alternative 
route to unravel manifestations of the intricate bound state problem in QCD. 

Among hadronic observables, generalized parton distributions (GPDs) [1] are the most elab- 
orate light-cone correlation functions containing simultaneous information on both position and 
momentum distributions of strongly interacting constituents [21 [3]. Similarly to conventional 
coUinear parton densities which are measured in deeply inelastic scattering, GPDs can be probed 
in experiments involving electroweak bosons with production of a real photon or mesons in the 
final state (see Ref. [1] for reviews). While the meson production receives contamination from 
additional hadronic final states, the former are free from these uncertainties and provide a clean 
access to hadronic inner content through GPDs. 

To date, the most complete analytic calculations were performed for photon leptoproduction 
cross section off spin-zero and spin-one-half targets and were limited to twist-three accuracy. 
This approximation implies that only terms suppressed by a single power of the hard photon 
virtuality were kept in all analytical expressions, corresponding both to leptonic and hadronic 
parts of amplitudes. The latter is defined by the expectation value of the chronological product 
of two electromagnetic currents between in and out hadronic states. 



This tensor is parameterized in terms of the so-called Compton form factors (CFFs) Q) 
which enter as coefficients in front of independent Lorentz structures. The CFFs depend on the 
generalized Bjorken-like scaling variable ^, the squared momentum transfer t, and the photon 
virtuality qf = —Q^. The QCD factorization theorems are indispensable in separating CFFs in 
terms of short-distance coefficient functions C(x, ^; controllable via conventional perturba- 

tion theory in QCD coupling constant, and long-distance dynamics encoded in GPDs F{x, ^, t; /i). 



The state-of-the-art considerations of the hadronic tensor (11.11) were done in the twist-three 
approximation El [TJ El E] • The hierarchy of hadronic matrix elements of higher twist operators 
emerging in the operator product expansion of Eq. (II. ip suggests smallness of their effect on event 
rates even at rather low virtualities. This phenomenon is well-known in deeply inelastic scattering 
and yields precocious scaling of corresponding observables making the neglect of operators of 
twist-four and higher legitimate. On the other hand, the approximation to merely leading and 
first subleading contributions stemming from the leptonic tensor were an artifact of matching the 
expansion of hadronic and leptonic parts. When the latter approximation is waived, numerical 
considerations demonstrate significant deviations between the two predictions for the kinematics 
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Figure 1: Amplitudes contributing to the photon leptoproduction cross section. The first one 
(i.e, the left-most) is the DVCS amplitude factorized into GPDs while the other two are the 
Bethe-Heitler amplitudes parameterized by hadronic electromagnetic form factors. 

of Jefferson Lab experiments. Therefore, in the present paper we perform a refined analysis 
starting with leptoproduction cross section of real photons off a spinless target. 

Our subsequent consideration is organized as follows. In the next section, we present a 
brief profile of the twist-three formalism we have developed in our earlier work [10]. In Section 
m we introduce the formalism of helicity amplitudes pTl [T2] which proves to be very efficient 
in separating power suppressed effects arising from strong-coupling dynamics in the form of 
higher-twist correlations, on the one hand, and kinematical effects due to nonvanishing masses 
of hadrons and momentum transfer in the t— channel, on the other. As previously, the hadronic 
part is calculated to twist-three accuracy, while we account for aforementioned kinematically- 
suppressed contributions exactly and thus keep all power effects in the leptonic part. As a 
consequence, the very transparent classification scheme of Ref . pjO], which allows one to identify 
Fourier harmonics in the azimuthal angle with specific twists of contributing GPDs, ceases its 
existence at small photon virtuality in the valence region. Then in Sections 14.31 and [5l we provide 
an estimate of different contributions to the cross section to get a handle on the most sizable 
effects. Finally, we conclude. The discussion of the kinematics is deferred to Appendix \^ while 
the explicit expressions for leptonic helicity amplitudes are summarized in Appendix [Bl 

2 Electroproduction cross section 

The main focus of our present analysis is the four- fold cross section for scattering of a light 
lepton £ = off a spinless hadron h and production of a photon in the final state, i{k)h{pi) —>■ 

e{k')h{p2Hq2), 

a^XBV 
da = , 

The phase space of the process is parameterized by the Bjorken variable xb = Q^/(2pi ■ gi) 
defined by the virtual photon euclidean mass = —qf of momentum qi = k — k', the squared 
momentum transfer t = A"^ with A = P2 — Pi, the lepton energy loss y = pi ■ qi/pi ■ k, and 
finally the azimuthal angle of the outgoing hadron. The dependence on the latter provides a 
very important handle on different combinations of twist-two and -three GPDs which enter the 
hadronic amplitudes fll.ll) . The cross section depends on small kinematical parameters which we 
will account for exactly in the present consideration. One of them has already appeared explicitly 
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dxBdyd\t\d(j) . 



(2.1) 



2 



in Eq. (I2.ip and is given by the ratio of the hadronic mass M to the photon virtuahty Q, 

M 

e = 2x^—. (2.2) 

The other will be introduced below. 

According to Fig. [T], the amplitude of the process T is a sum of two distinct contributions, one 
involving the deeply virtual Compton scattering (DVCS) tensor (11. ip and termed T^^^^ and the 
other one with leptonic Bethe-Heitler (BH) subprocess coupled to the hadronic electromagnetic 
current parameterized via the (pseudo) scalar form factor F{t) as 

J, = {P2\3M\Pi) = iPi+P2)^Fit) , (2.3) 

and dubbed T^^. The latter is real (to the lowest order in the QED fine structure constant) and 
F{t) is taken from other measurements. The azimuthal angular dependence of each of the three 
terms in 

^2 _ |^BH|2 _|_ |^DVCS|2 _|_ J {2 A) 

with the interference term 

arises from the Lorentz-invariant scalar products defining the leptonic and hadronic parts of 
amplitudes as explained at length below. Since the square of the Bethe-Heitler amplitude was 
computed exactly in Ref. p!3], we will focus our attention on the remaining contributions. Both 
of them are expressed by contractions of leptonic tensors with corresponding hadronic transition 
amplitudes, involving either the square of the DVCS amplitudes for |7"DvcS|2^ being linear 
both in DVCS and hadronic electromagnetic current for X, 

I^DVcsp ^ _^[^,.(X)T^^(j,r^y _ (2.6) 

^ = A'V^M {^-(^>^^''('''')* + t-l^rmm-r} . (2.7) 

Here the + (— ) sign in the interference stands for the negatively (positively) charged lepton 
beam. The leptonic tensors to the lowest order in the fine structure constant reacll], respectively, 

L^u{>^) = (k^j,kl + kyk'^ - k ■ k'g^^ + iXe^ukk') , (2.8) 

L,ur{\) = ^^^^^^^^^^tri(l-A75){7.(^-^)"S. + 7.(r+^)"S.^'}7M^, (2.9) 

where the lepton mass has been set to zero. The rescaled BH propagators 

ik^_ 2k.A _ {k-Ay _ t-2k.A 
^1 = Q2 - ^ + ' ^2 = , (2.10) 

emerge as contaminating sources of the azimuthal angle dependence which interfere with Fourier 
harmonics accompanying the generalized Compton form factors if expanded in inverse powers of 
the large photon virtuality Q^. Thus they will be treated exactly. 



^We adopt the conventions for Dirac matrices and Lorentz tensors from Itzykson and Zuber [14], e.g., gO^^a _ 
+1. We assume that the lepton hehcity is positive, i.e., A = +1 if the spin is ahgned with the direction of the 
lepton three-momentum. 
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Figure 2: The kinematics of the leptoproduction in the target rest frame. The ^-direction is 
chosen counter-along the three-momentum of the incoming virtual photon. The lepton three- 
momenta form the lepton scattering plane, while the recoiled proton and outgoing real photon 
define the hadronic scattering plane. In this reference system the azimuthal angle of the scattered 
lepton vanishes, while the azimuthal angle between the lepton plane and the recoiled proton 
momentum is 0. 

To continue our discussion we choose the target rest frame as shown in Fig. [21 The explicit 
components of particle's momenta in this frame are defined in Appendix |Al Guarded with these, 
one immediately computes all invariant products in terms of variables of the phase space of the 
process. For instance, one finds for the angular dependence of the BH propagators, 

where the 1 / Q-power suppressed kinematical if —factor, also showing up below in the Fourier 
expansion (13. ip - (13.21) . reads 

4xB(l-XB) + e2 t' \ ^2 12) 

with the plus sign taken for the square root in Eq. (12.111) and the variable t' standing for 

t = t ^min • (2-13) 

The variable K vanishes at the kinematical boundary t = tmiru determined by the minimal value 
of the momentum transfer in t— channel 




4xb(1 - xb) + 
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as well as at maximal value of the lepton energy loss ?/max = 2e ^(vT+~e^ — 1). 

3 BKM approximation 

In the frame we have chosen for our analysis, the contractions between the leptonic and hadronic 
tensor structures yield finite sums of Fourier harmonics, whose maximal frequencies are defined 
by the the rank-m of the leptonic tensor in the incoming lepton momentum k^. Note, however, 
that the polarization-dependent part of the leptonic tensors possesses one power of the four- vector 
fc^ less than in the unpolarized sector. As a consequence, the highest harmonic accompanied by 
the lepton helicity A will be sin([m — !]</)) rather than sin(m(/)), such that one finds for squared 
amplitudes 



±e 
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The generation of new harmonics in the azimuthal angle terminates at the twist-three level. The 
leading term (in inverse powers of the hard scale Q) in each Fourier coefficiently was computed in 
Ref. [13] and will be presented below for the sake of comparison with improved approximations 
computed in the next section. Here we will point out that the leading contributions to cf , sf as 
well as Cq^^^ emerge from twist-two GPDs [121 US]- The rest of the Fourier harmonics provides 
an additional angular dependence and is given in terms of twist-two, i.e., for cj, and twist-three 
GPDs, i.e., for cf^*"^, sf^*"^, cf, and sl^. The harmonics proportional to cos (30) [cos (20)] or 
sin (30) [sin (20)] in the interference [squared DVCS] term stem from the twist-two double helicity- 
fiip gluonic GPDs alone. They are not contaminated by any twist-two quark amplitudes, however, 
will be affected by twist-four power corrections ^16j . We neglect in our consequent considerations 
the effects of dynamical higher-twist (larger than three) contributions by the token alluded to in 
the introduction. 

3.1 Squared DVCS amplitude 

The Fourier coefficients of |7"Dvcs|2 naturally split into the product of factors depending on the 
leptonic kinematical variables and hadronic functions C^^*"^ 



^Dvcs _ 2(2 -2y + y'^)C^^^^ {H, H*] Ht, Ut) , (3.3) 

- Tn All ,9m V ' 



.DVCS 
-1 

.DVCS 



The latter are bilinear in the CFFs and respectively read 

CDVCS f n_/ n/^ nj nj^ \ njnj^ i nj 
[n, n , rtT, rLT) = '^'^ + 7773 u'^TrLj^ , 



^We adopt here the notation of Ref. pLO^ rather than of Ref. [13] 
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c^^^s (T^efi; n*; n*T) = (n* + jp^^^^T | (s.e) 

CDVCS /I-/ 1-/* \ 
2^ [rt, rtrp) = rLrtj- . 

At this point, we would like to recall that the spinless hadron acquires only two types of leading 
twist GPDs, unpolarized quark GPD H and gluon transversity GPD Ht giving corresponding 
names to the CFFs fll.2p . In the adopted approximation, the (^^^^ harmonic is expressed via the 
twist-two OFF 7i, while the coefficients cf and s^^*-^^ arise from the interference of twist-two 
and effective twist-three CFFs, 

= -2^ (yT^'^ + - 'H-^ , (3.7) 

with the CFFs related to functions given by a convolution of the twist-two GPD H 
and the so-called Wandzura-Wilczek kernel provided one neglects dynamical quark-gluon-quark 
correlation functions [13]. Here the generalized Bjorken variable ^ is related to the usual one xb 
via C, — xb/(2 — xr). Finally, the Fourier coefficients c^^^^ and s^^*"^ are induced by the gluon 
transversity CFFqfl. 



3.2 Interference of Bethe-Heitler and DVCS amplitudes 

Next, let us present approximate results for the interference term which is the most promising 
observable for the phenomenology of GPDs since it is linear in CFFs. This property simplifies 
the disentanglement of CFFs from experimental measurements. The Fourier harmonics have the 
form: 



( 2-y) 
l-y " ' Q 



4 = -8(2 - y)^e\ V^^' + 7=^^(1 - l/)(2 - ^b) (H) , (3.8) 





8AJ-(2-2y + y2)3?e- 

IQK^ [-(2-y) 3fJel 
2 — a^B 1 3^111 J 



(W) , (3.10) 



^eC^inr) ■ (3.11) 



M2(2-xb)2 



Explicit calculations demonstrate that the twist-three harmonics, i.e., cf and sf have the same 
functional dependence as the twist-two coefficients, cf and . However, this is not the case for 
cf and sf which emerge due to gluon helicity-ffip CFF Ht, 



c^{n) = Fn, c^ = FnT. (3.12) 

This set of formulas forms the complete result for the real-photon leptoproduction cross section 
in the twist-three approximation. 



•^These were omitted in our earlier consideration [13] of DVCS on spinless targets since they enter the amplitudes 
suppressed by a power of the QCD coupling constant. Note that Ht also enters as an {as/n)'^ suppressed 
contribution to cjj'^*^^ and as a as/ir suppressed effect in the twist-three harmonics cf^*-'^ and s^^*-^^, cf. Eq. 

(EH). 
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4 Helicity amplitudes 



The consideration of the previous section was restricted to the twist-three approximation for the 
dynamical as well as kinematical effects. Higher-twist operator product expansion analyses of the 
off-forward Compton scattering amplitude, similar to the one performed for the deeply inelastic 
scattering [17], is intrinsically involved due to complications and ambiguities in the choice of 
operator basis. The incorporation of kinematical power-suppressed effects is relatively straight- 
forward. In order to achieve this in the most efficient manner we separate power corrections that 
arise from the leptonic and hadronic parts by evaluating photon helicity amplitudes utilizing the 
polarization vectors for the incoming and outgoing photons in the target rest frame as defined 
in Appendix El in Eqs. (1A.6I) and (lA.Tj) . respectively. 



4.1 Squared DVCS term 



Using the completeness relations ( lA.lOl) and ( lA.llI) for the photon polarization vectors, we can 



rewrite the square of the DVCS amplitude (12. 6p as follows 

^ a=-,0,+ f)=-,0,+ 

in terms of helicity amplitudes labelled by the helicity states of the (initial) photon. These are 
defined by contractions of the Lorentz covariant amplitudes with the photon polarization vectors 
flX6|) and flAlTl) . 

Cab{X,^)=e't*ia)C,,{\)eUb), (4.2) 

and 

= i-ir-'er{c)T,.eUa) , (4.3) 

where the phase (— 1)"^"^ takes care of the signature factor in the completeness relation (lA.lOl) . 

The helicity amplitudes (14. 3 p are constrained by the parity conservation and, as a consequence, 
we have just three independent functions, 

-T-DVCS _ ^DVCS -T-DVCS _ -7-DVCS -7-DVCS _ -7-DVCS /a a\ 

^++ ~ ^ — ) ^0+ ~ -'O" ' ~ ■'+- ■ l^-^J 

The leptonic helicity amplitudes are calculated exactly from the definitions (12. 8p and (1A.6P and 
read 



with the remaining ones related to already found by parity and time-reversal invariance 
/:o-(A,0) = A+(-A, -0) , /:±,o(A,0) = >Co,±(-A,0) , 

/:__(A) = £++(-A) , /:_+(0) = £+-(-0) • ^ ' ' 
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Using the above relations among helicity amplitudes, we can cast the squared DVCS amplitude 
in the following form 

g2|^Dvcs|2 ^ 2/:++(A = 0) [T++ (r++)* + r+_ (T+_)*] + 2/:ooTo+ (ro+)* 

+ [£o+(A, 0) + Co+i-\, -0)] To+ [r++ + r+_]* 
+ [£o+(-A, 0) + Co+{X, -4>)] [T++ + T+_] (To+)* 

+ [/:+_(A,0)+£+_(-A,-0)][r++(T+_r+T+_(T++r] . (4.10) 

These findings immediately allow one to get the Fourier coefficients in the refined approximation. 
Substituting the hadronic helicity amplitudes computed to the twist-three accuracy from Lorentz- 
covariant parameterizations of the Compton tensor from Section 15.11 

T^^"^^ = n + 0{l/Q^), (4.11) 
To^vcs = -^^nf + 0{l/Q^) + 0{aJQ) , (4.12) 

n^=^-p^S«-+^'(vs^). (4.13) 

where the effective twist-three CFF Tif^ was defined earlier in Eq. (13.71) . we can read off the 
kinematically improved DVCS harmonics, 

Co =2 C (^>^.^^>^t)+(2-xb)^(1 + 62)' ^^-^^^ 

[ cfvcs ^ ^ f (2 - y) 3?e ) ^^vcs f^es. \ ^ ^^^^^^ 



Here the C functions are the same as appeared in Eqs. (13.61) . Thus, to restore the power suppressed 
contribution in the leptonic part it suffices to perform the following substitutions, found by 
comparing Eqs. ([S3D - dS^D with fimj) - fHTTB]) : 



2-2|/ + t/2 + 4^2 

^ (4.17) 



2-2y + y^ 



l + e2 

2-1/1 .1(2 



-Xy j 1 + e2 [ -XyVl+1 

g2^2 g2^2 



M2 M2(l + e2) 

4.2 Interference terms 

Finally let us turn to improving the interference X by treating it in the manner completely 
analogous to the previous consideration for the squared DVCS term. As a result, one finds 

^ = tvMV2{<pf ^^^ E {^ati\mt+{c^ti\mty} , (4.18) 
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where the leptonic part read 



/:f,(A,0) = eria)L,,riPi+P2relib) , (4.19) 

while the definition of the hadronic amphtude is self-obvious. In the present formalism, the 
angular dependence is entirely contained in the leptonic part and is rather intricate. Rewriting 
the interference in the form 

{C^^ + T++ + {C^^ + C^_) To+ + {C^^ + C^_) T_+ + c.c. , 

we introduce the decomposition of the leptonic amplitudes in Fourier harmonics encoding the 
azimuthal dependence, 

C^a^, + C^.a-b = -TT^ Ie C'„,(n) cos(n0) + a E Satin) sin(n0) | . (4.20) 

The explicit expressions for the exact Fourier coefficients in the leptonic tensor are given in Ap- 
pendix [Bl Notice that once one computes the leptonic part exactly rather than to the twist-three 
accuracy, as was done in Section [3l the simple one-to-one relation between Fourier coefficients 
and twist expansion in terms of CFFs is lost beyond the 1/ Q^— accuracy. This does not prevent 
one however to project out the real and imaginary parts of separate CFFs. 



±e^F{t) 



4.3 How robust is leading approximation for harmonics? 

As we will demonstrate in this section, differences between the approximate and exact amplitudes 
can be quite significant due to numerical enhancements of power-suppressed contributions in 
the valence and large— xb kinematical regions. We ignore the helicity-fiip amplitudes in our 
consideration and focus our attention on the consequences of the improvements in the twist-two 
sector encoded into the approximate Fourier coefficients and in Eq. (13.91) . At first let us 
propose a "hot fix" which is an approximation to the exact result and accounts for the most 
significant source of enhanced kinematical corrections, on the one hand, but leaves dynamical 
corrections out of the picture, on the other. This constitutes in the replacements 

f -(2 - 2j, + 1 f C±±(n = 1)1 

Xy(2-y) P1s,.(„=1)1 '"-'l' 

in Eq. (13. 9p . The admixture of higher harmonics proportional to Ti is not large for the present 
experiments, however, one should take care of the zero harmonics by the substitution in (13. 8p 

- 8(2 - y)|il^ir2 ^ _ ^)(2 - xb)| ^ C^±{n = 0) . (4.22) 

Notice that in this approximation the constant contribution, suppressed by the power of 1/ Q 
compared to the first harmonics, is entirely determined by the twist-two CFFs Ti. 

This "hot fix" provides a significant improvement of the leading approximation, as demon- 
strated in Fig. [3] where we compared it with exact formulas using for illustrative purposes the 
target with the mass M = 0.94 GeV . There, in the left and right panels, we present the plots for 
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Figure 3: Odd (a,b) and even (c,d) harmonics of the interference term (13.21) versus cj) as they arise 
from the hehcity non-flip distribution amphtude 7^+ (set to one) for JLab (a,c) and HERMES 
(c,d) kinematics as specified in Eqs. (I4.23P and (I4.24p . respectively. The dotted and dashed curves 
emerge from si sin(0) and cq + Ci cos(0) in the approximations of Eqs. (13. 8p - (13.90 and (I4.2ip - 
(I4.22p . respectively, while the solid ones contain all harmonics (I4.26p - (14.270 . 



the typical kinematical setups, with relatively large values of t', of the present JLab {Ef. = 5.7 
GeV), 

^' = -0.3GeV^ xb = 0.3, Q^ = 1.5GeV% (4.23) 

and HERMES {E^ = 27.5 GeV), 

^' = -0.3GeV^ xb = 0.1, Q'^ = 2.5GeV\ (4.24) 

experiments, respectively. A naked eye inspection of the plots of the low— Q/high— o^b JLab 
kinematics exhibits the evident feature that the approximation (13.80 and (13.90 (shown by dotted 
curves) provides a sizeable false enhancement effect of the leading twist harmonics compared with 
the exact result (solid lines). Fortunately, the contamination by higher harmonics is small and 
so the substitutions (14.220 - (14.210 provide already a good agreement with exact results (dashed 
line in Fig. [3]). For low— Q^/large— xb, we find ~ 30% and ~ 70% deviations for the first odd and 
even harmonics in the angle 0, respectively. For HERMES, where xb is smaller and larger, 
the approximations (13.80 and (13. 9p are justified as it is obvious from Fig. [31 We note, however, 
that also for this kinematics at larger values of xb one may find relative deviations from the exact 
results of the order of 15% or so. 

One should also be concerned about the size of twist-three effects, originating from the 
longitudinal-transverse helicity amplitude Tq±. They are formally suppressed by ^J—t'/Q and 
show up in higher harmonics in the leading approximation of Section [3l which, however, are 
kinematically contaminated by lower ones when computed exactly. For instance, for the JLab 
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kinematics ( 14.23^ we find 



J oc - 1.34 [cos(0) + 0.51 - 0.14 cos(2(/))] ^eH - 0.68 [cos(20) - 0.32 + 0.21 cos(0)] 3fJe Uf 
+ 1.16A[sin(0) + 0.04 sin(2</))] 3m 7^ + O.23A[sin(20) + 0.79 sin(0)]$5m7^f. (4.25) 

If tfie twist-tree CFF Tif^ is comparable in magnitude to tlie twist-two Ti, it may provide a 
sizeable contribution to the lower harmonics, i.e., sin(</)), cos(0) and the constant part. The size 
of the sin(20) and cos(20) harmonics will serve as an accurate estimate of this admixture. Note 



that for (very) small values of ^—t'/Q, the twist-two and -three harmonics start to decouple. 
However, even for HERMES kinematics one must account for these mixing effects for average 
values of -t ~ 0.3 GeV^. 

Finally, we mention that in the present estimate (14.251) we ignored the transverse-transverse 
helicity-fiip amplitude T+_ since it safely decouples from the other ones. It shows up in the 
cos(30) harmonics, which is only slightly affected by the twist-two and -three CFFs. On the 
other hand, its contribution to the lower harmonics is negligible. 

Since the ultimate goal of the exclusive cross section (12.11) measurements is the extraction 
of CFFs and consequently clean separation of GPDs, which shed light on the internal structure 
of hadrons, data analyses have to disentangle the intricate dependence of kinematical factors 
on dynamical variables from the one of GPDs themselves. Therefore, the role of 1/Q power 
suppressed contributions should not be underestimated as these enter dressed with numerically 
enhanced factors being function of y, xb and t. Therefore, for the most robust CFF extraction, 
we suggest to utilize the following general formulas for the Fourier coefficients in data analyses 



cl = C++(n) JJeC^ (H) + | Co+(n) 5JeC^ (uf) + ^^^^ ^ C.^in) ^eC^ {Ut) 



(4.26) 



si = S^^in) 53m (H) + ^ SoM ^SmC^ (nf) + ^ ^ 5_+(n) 53m (Ht) , 

in the interference term (13.21) . Though this results into a more involved extraction procedure, 
compared to that one described in Ref. [lOj, it will not introduce any complications of principle. 



5 Uncertainties arising from the twist expansion 

There are various possibilities to write down the tensor decomposition for the DVCS tensor (11.11) 
respecting its gauge invariance and Lorentz covariance. Each of them will be associated with a 
specific set of CFFs. Of course, all parameterizations are equivalent. However, relating physical 
CFFs to partonic GPDs requires an expansion in inverse powers of the hard scale Q, known as 
the twist expansion. As a consequence, expressions obtained for different parameterizations will 
only be equivalent to the considered order in the 1/ Q— expansion and differ when higher order 
terms are taken into account. Let us address several Lorentz structures used in the literature 
and discuss their manifestation in physical observables. 

Presently, the hadronic Compton amplitude (II. ip has been worked out to the twist-three 
accuracy within the GPD formalism El El El E] • While the leading twist sector is known in the 
next-to-leading [18] and next-to-next-to-leading [19] orders in the modified minimal subtraction 
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and the special conformal subtraction schemes, respectively, the twist-three coefficient functions 
are available in the handbag approximation only (with a partial one-loop calculation reported 
in Ref. |20|). The restriction of the current analyses to account for 1/Q suppressed effects is a 
drawback which yields insufficient accuracy to render the DVCS tensor gauge invariant and obey 
the transversality condition, arising from the current conservation, 

= ^2%, = 0, (5.1) 

exactly. To illustrate the interplay of higher twist contributions and the restoration of the 
transversality condition of the DVCS tensor, let us recall the role of the twist-three effects in this 
endeavor. For a scalar target, the twist-two result in the handbag approximation reads: 

T^^) = ~~r^ ■ (l)9fiu - q^iPu - quP^, - iP^Pu] T-t^i, t; Q^) , (5.2) 

where we used the symmetric variables p = pi+p2 and q = \{qi + q2) which spawn the generalized 
Bjorken variable ^ = Q'^/p ■ q defined in terms of the average photon virtuality = —q^. In the 
non- forward kinematics of DVCS, obviously this tensor does not respect the current conservation 
and so one finds after contraction, e.g., with qi^, that the transversality condition (15. ip is violated 
by A;^ = Ay—r]py which, according to the power counting, is a twist-three effect. The variable r] is 
the t-channel scaling variable known as skewness t] = {A-q)/{p-q). It is related to the generalized 
Bjorken variable ^ for the DVCS kinematics via the relation r] = — ^/(H-t/2Q^) ^ — ^. Evaluation 
of the twist-three amplitudes yields beside a contribution proportional to the function H3 

Tj^J = ^ [Ai (g2. + ^P.) + (gi^ + ^p,)] , (5.3) 

also an antisymmetric tensor accompanied by the twist-two CFF Ti 



1 



H 



p-q 

The latter when combined with Eq. (15. 2p . T^^^ +5T^'^\ generates the replacement of the twist-two 
Lorentz structure with 

-{P ■ <l)9t^u + q^lPv + quP^Ji + ^P^,Pv ^ -{p ■ q)9^lv + ql^,Pv + q2vP^, + ^P^,Pu ■ 

It is straightforward to check now that the violation of the transversality condition is postponed 
to the twist-four level, i.e., C(Q~^). Thus addressing it requires a full-fledged twist-four analysis 
along the lines of Ref. [1^. Since such a calculation has not been performed yet in a consistent 
manner, we will take a pragmatic point of view and restore the exact gauge invariance in several 
superficially inequivalent ways and study the resulting numerical differences. 

5.1 Lorentz decomposition of the Compton tensor 

The complete twist-three result for the DVCS tensor for a scalar target has three independent 
CFFs, which is the maximal possible number allowed by the underlying symmetries. The result 
reads in terms of physical four-momenta as follows 

1 / A^ \ 

+ -LA^^A^,7^T + [Ax^ {q2u + iPu) + A^, (gi^ + ip^)] . 
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As mentioned above the two Lorentz tensors, proportional to the CFFs Ti and Tis, are transverse 
up to twist-four terms. The third structure is a symmetric and traceless tensor, which is induced 
by the gluon transversity (photon hehcity flip by two units) [2T1 ES]- It is only transversal to 
the order A_|_/M. 

Utilizing the equation of motion at twist-four level, should allow one to satisfy the transver- 
sality condition ( 15.ip exactly, since higher than geometric twist-four contribution are formally 
absent in the handbag diagram being proportional to derivatives of Dirac delta function. The 
restoration of the transversality condition for the Lorentz structures, proportional to Ti and 
7-^3 requires a t— dependent term, since neglecting such terms results in a failure of (15. ip [21]. 
The "minimalist" restoration in the twist-two sector requires to add a {t/Q'^)p^Pu contribution, 
leading to the replacement 



(p ■ Q)9f,u - qif,Pu - q2uPf, - iPixPu 



(p • <l)9, 



4Q2 



^Pi^Pu- (5.4) 



In the twist-three sector one needs a 

t q^Pv 



p-q 



proportional term, which can be absorbed by adding power suppressed contributions to the 
following four-vectors: 



A 



2Q' 



At ) 



t 



4Q2 



(5.5) 



with ± sign in the first equation standing for indices associated with outgoing/incoming photons. 
The restoration of the gauge invariance for the gluon transversity is more cumbersome, since it 
will already affected at the twist-three level. It can be achieved by adding a proportional 
terms [16]. For later illustration we adopt here the recipes (15.41) and (15.51) . Finally, the DVCS 
tensor can then be written as 



T, 



1 



+ 



p-q 

1 

1 



(P ■ q)9^,u - qif,Pu - q2uP^ - U + 



4Q2 



^P^iPu 



(5.6) 



p-q 



2Q^ 
q2v + ( 1 + 



2Q 



4g2 



^p. Ax. + 



2Q' 



;^Pu 



4Q2 



^Pn 



A,„ - 



2Q- 



;iPv 



For DVCS kinematics this parameterization is complete and, thus, it can be uniquely mapped 
into a different form, e.g., used in Ref. [25] . 

Another way to restore gauge invariance of the amplitude (11.11) beyond twist-three accuracy 
is by introducing a projector [7] 



q^ 

qi - q2 



(5.7) 
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fulfilling the conditions V^^qiy = q2fiV^^ = 0. They provide a transverse Compton tensor when 
contracted on both sides with the twist-two DVCS amplitude Tj^J . When expanded to twist-three 
accuracy this procedure reproduces gauge-restoring terms coinciding with the ones obtained from 
the explicit analysis reviewed in the preamble to Section [51 Notice however that contrary to the 
consideration in the previous paragraph, this recipe generates an infinite tower of higher-twist 
contributions when expanded in inverse powers of the average photon virtuality. 

Finally, let us address the formalism of Ref. [26]. In applications of the QCD- improved 
parton model to high-energy scattering, it is customary to parameterize the Compton tensor 
in terms of the light-cone vector and its tangent such that = = and n ■ h = 1. 
In particular, these can be chosen as the plus and minus components of the initial photon and 
proton momenta, respectively, cf. Fig. O However, the choice of the light-cone vectors is not 
unique, since it implicitly refers to a given reference frame. For instance, the parameterization 
used in Ref. I26l. 



4M2 



t 



-n 



- 4M^ — t 

= -(20n^ + e z + A±a* , (5.8) 



is done in a "collinear" frame {pi± = —p2±)- We note that the scaling variables ^' and ^ are 
proportional to the generalized Bjorken variable ^ and skewness t] ~ — ^, respectively, but, differ 
from them by power suppressed corrections. The DVCS tensor parameterized in terms of the 
light-cone vectors rather than physical four-vectors reads [26] 



T, 



+ 
+ 



n ■ q2 



H 



Ax^ + Al 






h ■ q2 




A^^A^, 


Ai 


( 


M2 


2M2 I 



%A_Li. 

n ■ q2 



(5.9) 



n 



T • 



In the first and third Lorentz structure the linear term A±y exactly restores the transversal- 
ity, while for the second Lorentz structure a A^— proportional term is needed. Note that the 
restoration of transversality for the gluon transversity-induced Lorentz structure essentially dif- 
fers from the prescription employed in the parameterization fl5.6p . In the frame we are adopting, 
the transversal helicity amplitudes, 7^ 7^, are given by CFFs Ti and Tix, while the the 
longitudinal one, 7^ — >■ 7t, is entirely related to the twist-three CFF Hs. 

To get rid of the frame dependence in the parameterization (15.91) . we use the relations (15.81) 
to express the DVCS tensor in terms of the physical momenta. After this, the tensor takes a 
frame-independent Lorentz covariant form. Comparing it to the parameterization (15.61) we can 
read off a rather cumbersome relation among the two sets of CFFs. 



5.2 Numerical estimates 

The leading contribution to each hadronic helicity amplitude from any tensor decomposition, 
either (15.41) or (15.61) or (15.91) . is universal and is given by Eqs. (14. lip - (14.130 . Differences will 
arise starting with l/Q^ contributions. To illustrate uncertainties related to the twist-four effects 
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we evaluate the helicity amplitudes to the order The differences in the helicity conserved 

amplitude from the calculated (15.41) and improved (15.61) DVCS tensor is 



^++ lEq.lEl " ^++ lEq.lESl - ^ ( I • 

It can be considered as negligible and plays practically no role. The ambiguities in the restoration 
of transversality with (15.41) shows up mainly in the gluon transversity induced sector yielding 

^Dvcs, „ ^Dvcs, „ (4M^ - t)xi 4(1 - xb)xbM^ - (4 - 3xB)t ^ ^n(Mn^\ 

^++ lEq.|53}-^++ lEq.l53} = Q2 (2-Xb)^M2 ^T + C(l/Q). 

(5.11) 

It is suppressed for small xb by a factor x%. The CFFs and 7^3 enter here as a x^tt'/ and 
x^M'^t' / suppressed contributions, which are practically very small. 

For the longitudinal-transverse helicity amplitude, the higher twist effects are more pro- 
nounced and yield the difference 

T-DVCSi , , T-DVCSi , , fKTO\ 



xbKI 



n 



2v^(2-xb)Q3 V L ^ / 2-xb 



1 -2x 



1 + 0(1/Q') -2 ; "•^^ 7^3 1 + 0(1/Q 



ir 2M24-t(2-2a;B + 4) .^ 
Q \/2M2(2 - xb)3 ^ 



1 + 0(1/Q 



2^ 



While the two recipes (El]) and ^M) yields xb(4M2 - t)7^f /4Q^ x|(4M2 - t)7^'=f^/4Q2 and 
a;|(4M2 — t)TCT/4:Q suppressed differences. 

Finally, for the transverse-transverse helicity-fiip amplitude, we find that the ambiguities due 
to kinematical higher- twist corrections enter at order l/Q^ for the recipe (15. 6p . i.e., 

r^'Wm - ^-:'\.m - o{i/Q') , (5.13) 

however, only at order 1/Q^ for the prescription (15.91) . 

To give a numerical example demonstrating the contamination of the leading contribution 
by the ambiguities of the power suppressed effects, let us start with the JLab kinematics (I4.23p . 
Numerically, we find that the non-flip helicity amplitude can be safely approximated by the CFF 
n, cf. Eq. fimi) : 

^Dvcs _ ) Q ggg Kn+i 0.011 \'H3 + i 0.019 ^Ht for <' Eq. (EH) \ . (5.14) 
[ 1.003 J [ 0.008 J [ 0.000 J [ Eq. (E3]) J 

Certainly, here we can practically set Tj^^'^^ = Ti. Unfortunately, for the longitudinal-transverse 
helicity flip amplitude we find rather strong deviations from the leading approximation (14.121) 
that are caused by kinematical corrections: 

(2-xb)QToD^^s 



r 0.997 ^ 




' 0.010 ^ 




' 0.019 ^ 




' Eq. (15. 41) ] 


0.996 


H 


0.011 




0.019 


>nT for < 


Eq. (15JD 


[ 1.003 ^ 




^ 0.008 ^ 




^ 0.000 ^ 




[ Eq. (15.9!) J 











' 0.16 ^ 




r Eq. (i5:iD ] 








hi 


0.34 


>nT for < 


Eq. ^ 


I 0.91 ^ 




0.02 ^ 




. 0.03 ^ 




^ Eq. (ED J 



(5.15) 
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Although the effect of restoration of the transversahty from the twist-two amphtude (15.40 is 
small within the recipe (15.61) (except for the gluon transversity) , the numerical values deviate 
considerably for (15. 9p . We suggest to rely for simplicity on the leading approximation (I4.12p . 
Finally, for the transverse-transverse helicity flip amplitude we again observe that the deviations 
from the leading approximation (14.131) are negligible except for the parameterization (15. 9p : 

(2 - X )2MTDVCS f 1-01 ] r -0.02 ] ( 0.06 ] f Eq. (El 

^ -+ ' ^ ^ T^T + <! -0.02 }n-{ 0.07 } H3 for { Eq. (TO 

(5.16) 



[ 0.81 J [ 0.02 J [ 0.04 J [ Eq. 



It is rather obvious that for decreasing xb and/or increasing the "kinematical" power 
corrections are getting reduced. For instance, for HERMES kinematics (14.240 we find for the 
most problematic longitudinal-transverse helicity-fiip amplitude 



(2-a:B)Qro°vcs 



V2K 







' -0.02 ^ 




' 0.05 ^ 




' Eq. (I5.4D ] 
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>nf + < 


-0.01 




-0.04 


> Ht for < 


Eq. f|5.6|) 






0.00 ^ 




. -0.02 . 




Eq. (15JJ) J 



(5.17) 

This exhibits the legitimacy of the leading approximation (14.121) employed in Refs. [TOt [13]. 



6 Conclusions 

The main goal of the present consideration was understanding of the power-suppressed effects in 
DVCS observables stemming from the exact account for kinematical contributions in the hadronic 
mass and momentum transfer t. Using the photon helicity amplitudes, we separated the 
leptopro duct ion cross section in terms of the leptonic and hadronic helicity amplitudes. The 
choice of the target rest frame with the z— axis directed (counter) along to the virtual photon 
three-momentum allowed one to localize its dependence on the azimuthal angle to the leptonic 
part. These were then computed exactly to leading order in QED fine structure constant thus 
improving approximate results of previous considerations p!0| fT3]. 

Numerical estimates performed for the current kinematics of JLab experiments demonstrated 
that, due to rather low virtuality of the hard photon and valence-region values for the Bjorken 
variable, the restriction to merely the leading approximation of Ref. fiUi [13] yields significant 
overestimate of event rates compared to the exact treatment. However, for higher values of the 
hard scale, typical for the HERMES experiment, the approximation of the earlier work becomes 
legitimate. We proposed a set of formulas for refined analysis of DVCS observables. Although, in 
the improved approximation the classification scheme of Ref. [10], according to which the Fourier 
harmonics are strictly associated with the twist of the contributing CFFs, is altered, this does 
not represent a difficulty of principle to extract CFFs from experimental observables. However, 
obviously the inversion problem becomes more tedious. 

Let us point out that the choice of Lorentz invariant kinematical variables in the evaluation 
of CFFs from the corresponding GPDs is also not unique as it has a cross talk with higher- 
twist contributions. The optimal choice should minimalize 1/ Q^— suppressed contributions. This 
problem was not the focus of our present analysis, where we used as in our previous studies a 
legitimate choice ^ = —rj = xb/(2 — xb) and took the photon virtuality as the large scale. 
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The consideration of the present work can be generahzed in a straightforward fashion to 
targets possessing nonvanishing spin, with nucleon being the most interesting one. We anticipate 
that our current analysis will quantitatively hold for DVCS off an unpolarized proton target 
too. The "hot fixes" fl4.17p . (14.211) . and (14.221) can be immediately used to improve on Eqs. 
(43) - (45), (54), and (53), respectively, of Ref. (TO]. Note that for the spin-one-half target new 
combinations of CFFs will emerge, e.g., AC in the interference term. However, we expect that 
such contributions, induced by the helicity flip of the outgoing proton, are relatively unimportant 
for the kinematics of present experiments. We will report on this analysis in future work. 

We are indebted to H. Avakian, M. Gargon, M. Guidal, and F. Sabatie for discussions which 
initiated our studies. This work was supported by the U.S. National Science Foundation under 
grants no. PHY-0456520 and no. PHY-0757394 and funds provided by the ASU College of Liberal 
Arts and Sciences. D.M. would like to thank the Particle Physics and Astrophysics Group at 
ASU for hospitality extended to him during the flnal stage of the work. 



A Kinematics in the target rest frame 

We flx our kinematics by going to the target rest frame and choosing the ^-component of the 
virtual photon momentum negative and the x-component of the incoming lepton being positive. 
The components of the corresponding four-vectors read 

pi = (M, 0,0,0), gi = -(l,0,0,-Vl + e2) , A; = — (1, sin^,, 0, cos^/) , (A.l) 

p \ /IIP 



2/e 



with the lepton scattering angle being 



,2 



1 + ^ t\ l-v - ^ 

cosei = -^^^, smei = ^-^J=^. (A.2) 

The outgoing momenta are parameterized in terms of the scattering angles in the hadronic plane, 
see Fig. [21 

Q2 = ~~Y7 — — (1, cos sin 6'-y, sin sin6'^, cos6'-y) , (A. 3) 

2iw 



»2 = [M ^, W-t + cos(/)sin6', J-t + sin(/)sin6', J-t + cos6' 1 , (A.4) 

^ I 2M'V 4M2 ^ 'V 4M2 ^ 'V 4M2 ' ' 

where the polar angles read in terms of the kinematical variables of the phase space 



The azimuthal angle of the photon ip^ is related to the one of the outoing hadron (p via (p^ = </)+7r. 

The explicit component form of the photon four-momenta allows one to construct their po- 
larization vectors: 

en±) = ^(0, 1, ±z, 0) , em = ^j^={-^lT7\ 0, 0, 1) , (A.6) 
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£0 ±) = ^ 0, , COS0 ± 2sin0, =F2cos0 H , sm0, — 

(A.7) 

which are defined up to an overall phase factor. The kinematical factor entering the last compo- 
nent of q2 reads 



and is related in an obvious manner to K of Eq. fl2.12p via 



/ e2 K 

K = \l^-y + jy^Q- (A.9) 

The photon polarization vectors obey the following completeness relations 

E e^,{h)er{h) - emeUO) = + ^ , (A.IO) 
h=-,+ ^1 



which are used in the main text to reduce the cross section to the product of helicity amplitudes. 



B Fourier harmonics in leptonic tensor 

Let us present explicit expressions for the Fourier coefficients entering the leptonic part of the 
interference term (14.181) . For the transverse-transverse harmonics we found 



- (TTTy [q^TTt^ ^ ' 



t 



+ 7^ (2-a;B) 1 + 



2 Wo ™ N 1 , I 2 ' 2xbJ Q 



^2- 



^ ) \ (2-XB)(i + yrT^) 



2 



C++(n = 2) 



1 - Vl + e2 + 3e^ XBt ' 

^ 1 + v/rT^-~ Q^, 



:i + e2)2 1 i + e2 + v/rT^Q2 



XBtt' f VTT^-i 

+ 1 - XB ^ + 



V 2 2x1 





= 3) = 


S++{n 


= 1) = 


S++{n 


= 2) 




4ri 
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1 + e2 [ 1 + e2 Q2 J ' 

(1 + e2)3/2 + VI + 6 2xs) I ^ ^ ^/^^ _ 2v/r+i^ /' 

while we got for the longitudinal-transverse ones, 



^o-^(n = 0) = ^^-^^ |e + 3 -| , (B.2) 



i2 



8^2(2 - -y-'^y^ 

(1 + e2)2 Q2 ' 

2 



So+{n = 1) 



^0 (n = 2) = ^ ^ 1 + - W 1 + , ^ 

Finally, the hehcity-flip transverse-transverse coefficients are 

r ( m 8(2 -y) f.o .Rox 
= °^ ^ (IT^^r " 2(1 + ^^-'^ 

+ ^r+i^ l'-"^ — ^ — +^j^r 

8^ f/o ^2 2- v^rT?/^v/rTe2-l + eV, M ^b*' 



C_,(„ = 2) = 4(2 - ,) (l - , - ly^ - 

J. A o , 2(1 - 3:b) 'i x^i' , , VI + + .Xb + (1 - J^b)^ i 



2 



C_+(n = 3) = -8X 1 - y - -y^ . .../2 ' 1 + 
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S-+{n 




a/TT?- 1\ XBt 

2^ ) 



I 



2xb) Q^J 
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